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Abstract. A discretisation of differential geometry using the Whitney forms of algebraic 
, topology is consistently extended via the introduction of a pairing on the space of chains. 

This pairing of chains enables us to give a definition of the discrete interior product and thus 
provides a solution to a notorious puzzle in discretisation techniques. Further prescriptions 
are made to introduce metric data, as a discrete substitute for the continuum vielbein, or 
Cartan formulation. The original topological data of the de Rham complex is then recovered 
■ as a discrete version of the Pontryagin class, a sketch of a few examples of the technique is 

also provided. A map of discrete differential geometry into the non-commutative geometry 
of graphs is constructed which shows in a precise way the difference between them. 

: 

In a wide range of problems in physics and engineering, computational challenges have become 
a major hurdle, for instance in fluid dynamics |Tj, electromagnetism |2j, and in fundamental 
i | physics. 

Many of the issues may be formulated as algebraic problems: how can we impose the prop- 
_C . erties of the continuum theory in the algorithm, in the definition of fields, or impose constraints 

J> \ on the evolution of the system, if we have in mind a symmetry of the system such as rota- 

k^J ' tional symmetry. As is known, the most concise way to formulate topological and geometrical 

properties is by use of the coordinate free language of differential forms to represent the fields 
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and use the mathematical operations on them. The latter are namely the wedge product, the 
exterior derivative, the Hodge star and finally the interior product (or contraction operation), 
respectively denoted by the symbols (A, d, *, ix) where X stands for a vector field. With this 
collection of operations we can define all operators of interest, such as the Lie and covariant 
derivatives. 

To define such rules on a lattice, or more precisely on a simplicial or hypercubic complex is 
the problem we address in this article. Much work has been done in this area, for instance by 
discretisation of the underlying topological theory [31 @J [51 El H[ and introducing relevant 
operators such as the discrete analogue of the contraction operator [21 E3 HTj . even though 
some algebraic properties are lost in the process. This structure is often referred to as a discrete 
exterior calculus. Of course, which properties are relevant will depend on the application at hand, 
for example in the electromagnetic theory, the duality between the electric and the magnetic 
fields is important. To find the optimal discretisation from an algebraic point of view is the 
problem of interest here. 

We will develop further our earlier proposal . In Section [21 we define an extension of 
the Whitney map, which maps discrete p-dimensional objects (p-chains) to p-forms, to W, 
which are maps from the Cartesian product of chains to forms, in such a way that the interior 
product is well-defined and exact . A key construction is to use the map W introduced to 



1 in the sense that from the continuum viewpoint, the result of the operation is correct. Not to be confused 
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give an alternate definition of two of the basic operations of differential geometry, namely the 
discrete wedge and Hodge star. In Section |SJ we will propose a new approach to discretising 
the metric in the same spirit although additional prescriptions are required, in Section a 
map between our discrete differential geometry and Non-commutative geometry of graphs is 
constructed this clearly shows the different between these two structures. The examples provided 
here are for two dimensions. In Section 5, we discuss the applications we are currently pursuing, 
namely, the Reidemeister torsion calculation in topological field theory, 2D lattice gravity and 
the discretisation of Maxwell's equations. 



with exact forms of algebraic topology 
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1 Introduction 

1.1 Importance of the Homology ring 

First, we recall some basic concepts of topology. The space of chains 

C*(K;R) = Q)C {p) (K;W) (1) 

p 

contain linear combinations of simplices with coefficients in the field R, on which one defines a 
collection of mappings, namely the co-chain space C*(K; R) which in turn, within the de Rham 
co-homology setting, may be identified with a finite subspace of the space of differential forms, 
called Whitney forms |13j (see the map W below). 

To our knowledge, all approaches currently considered to the discrete exterior calculus take 
all mappings from chains onto chains. To be more precise, the operations are all surjective in 
the space of chains, e.g, the discrete wedge product is defined as 

A K : C (p) (K; R) x C (q) (K; R) — C (p+q) (K; R). (2) 

The surjectivity in the definition of the discrete maps, we will think of as an assumption referred 
to as A. One can also appeal to the simplicial approximation theorem to motivate the latter. 

Because they were based on homological operations, for example, the dual complex L is 
constructed in order to define the discrete Hodge star (there are two Hodge star maps * K and 
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-k L mapping the complex K to its dual L and back), the cup product appears in the definition 
of the discrete wedge product; such mappings respect the homology ring 2 . 



Let CTup) denote the i-th p-simplex, and \f p ) coefficients in R. With the homology class of 
the chain 




written as [c(p)], then 

[°"(p)] ^ K W(q)] = [ a (p+q)] 



(6) 



Similarly, using the discrete Hodge star from K onto L 



,KrK] r L 1 
* [%)J = F(n-p)J' 



(7) 



where n is the dimension of the complex. The fact that the homology ring is captured does not 
however mean that the operations themselves are exact, not even that they satisfy their defining 
algebraic relations. We will see in the case of the discrete wedge that it is non-associative but 
in a way that is mild. 

The virtue of respecting the homology ring is easily understood, the operations give a chain 
that is in the correct homology class. That an operation is approximate does not alter this 
feature. Thus, enforcing this constraint from topology in the construction (which we will call 
the "Homological constraint", C), leads to a lack of exactness of the operations, apart from the 
discrete exterior derivative. 

The property of "exactness" denoted P below means that for a given simplex, the operation 
is numerically exact for any simplex as can be checked by application of the map W or W as 
we will see shortly. 

1.2 Lack of exactness of the operations 

In order to understand how this approximation appears concretely for the sake of the discussion, 
we enforce both the assumption A and the constraint C. To a given chain &i : (p), we associate a 
differential form W{ai^). The list of such forms in the case of the hyper-cubic complex can be 
found in Example 12.11 below. 

The immediate problem is that of exactness: take a Whitney form, say W([01]), while com- 
puting the exterior derivative diy([01]) gives a linear combination of Whitney forms and the 
exact result, there is no way to obtain closure of the Hodge star on the space of Whitney forms, 
i.e the list of Whitney forms does not include the Hodge duals of Whitney forms. In general, 
given a chain a, there is no chain a' such that -kW{a) = W(a'). Similarly, a discrete wedge, also 
to be defined below can only be approximate. From this simple observation, we conclude that 
a discrete wedge product or Hodge star satisfying A cannot be exact. 

Nevertheless, such discretisation has remarkable algebraic properties, since the discrete ex- 
terior derivative d K , which in the simplicial setting is the co-boundary operator denoted d K 
(which is motivated by the Stokes theorem Eq. I28J) is an exact operation and we do have 



dW = Wd 



(8) 



2 The cup product as a discrete wedge, and the discrete Hodge star satisfy respectively 



H( P )(K; R) U H( q )(K\ R) S H (p+q) (K; R), 
H{ P )(K\ R) S H {n _ p) (K;R). 



(3) 
(4) 
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Although the discrete wedge product is numerically approximate, the Leibniz rule is satisfied, 
and similarly, although the discrete Hodge star operator is approximate, it satisfies the correct 
algebraic properties and gives an approximate adjoint 

S K = (_i)«(p+i)+i d L * K (9) 

for the discrete operator d , which satisfies 

{5 K f = 0, (10) 

which may also be seen as a consequence of the exactness of d . The approximation coming from 
the other operation used, namely the discrete Hodge star is done within the correct co- homology 
class. 

To summarize, it has been found that under the assumption A, a collection of discrete maps 
for A , d K and * K and similarly on the dual complex L satisfying the constraint C is well- 
defined. Let us now turn to the question of discretising the interior product %x- 



1.3 Importance of exactness for the discrete interior product 

When one comes to the problem of supplementing the model with an interior product (the 
contraction operation), it is evident that: 

• The operation, defined from chains onto chains (i.e satisfying A) can only be approximate 
(violating P). 

• The algebraic properties do not hold, notably the anti-derivation property. 

The situation is similar to that of the discrete wedge product in that an approximation must be 
made, but in the present case, a discrete operator for i x using the operations (A , d K ,* ) will 
not satisfy exactness as was the case for d K , thus some algebraic properties are spoiled. Such 
proposal was initiated by Hirani jlUl II 1 j , who discretised the following continuum expression for 
the contraction operation on a n-dimensional manifold 

i xU W = (-\) k (n-k) + Jk) A ^ ( U ) 

where X b is the dual form to the vector field X and cjw is a fc-form. And so one may use this 
formula directly in the discrete setting with the discrete version of the operations A, *, d in place 
of the continuum ones. 

While the convergence of the operation in the continuum limit has been shown, from an 
algebraic point of view it is incomplete, due to the non-associativity of the wedge product 
alluded to above. One motivation for building such operator with more algebraic properties is 
to consider the Lie derivative given by the Cartan formula, 

L x = i x d + di x . (12) 

After defining the usual operations which allow to discuss topology in Section El we move 
on to the central construction, first, we introduce a larger space using a product on chains, this 
requires extending the Whitney map to this larger domain, the extension is not trivial and we are 
lead to check that the co-homology still holds. The natural one-to-one correspondence between 
chains and co-chains still holds while a larger subspace of differential forms is constructed. After 
the introduction of the interior product ix, we discuss its algebraic properties which are identical 
to the continuum ones, and derive the Lie derivative, arguing in what sense the Jacobi identities 
are piecewise satisfied. Then, following the logic of the introduction of the product space, we 
suggest a new Hodge star (denoted it E ) which does not require the introduction of a dual 



6 



V. de Beauce and S. Sen 



complex (this map was introduced in the discussion of lattice QCD in [20]). In the latter part, 
we present the simplicial version of the above, for completeness. 

In Section |31 we introduce metric in the present scheme, for this purpose, we find that the 
map W is essential. We introduce a discrete version of the vielbein and of the Cartan structure 
equations and then consider the Pontryagin class which leads to topological invariants. In 
Section HJ we discuss the relation of the present approach with the non-commutative geometry 
on finite sets as a discrete differential geometry proposed by Woronowicz (see also |24| ) . 
A relation is established after removing the smoothness associated to Whitney forms and to 
their generalization given here. We show how to extract the non-commutative setting from the 
present formulation, the two being not equivalent. A discussion of some applications to physics 
follows in Section 5. 

Before starting the construction, we ask the following question: what features of a manifold 
can we expect to capture on a simplicial complex ? 



1.4 Local versus global, simplices and manifolds 

Consider the triangulation of a n-dim manifold M by a simplicial complex K, in the sense of 
singular homology it can be identified as a combinatorial manifold and as a piecewise linear (PL) 
manifold, that is we have a manifold structure. In the standard simplex 

A n = {(x ,...,x n+1 ) : 0<Xi<l, Js j = l}cl n+1 l (13) 

i 

associated to each simplex, the continuous mappings a : A n — ► K are the singular simplices 
of K. Coordinates are thus provided, and it is in these coordinates that we will define the basis 
of Whitney forms and the extension defined here under the map W . 

Given two open sets U a and Up, associated to two neighbouring top degree simplices cr a r n ) 
and <Jp/ n )i the transition function 

<j>a/3 = <j>a O (pp 1 (14) 

defined on U a D Up is a homeomorphism. However it is not a diffeomorphism and so the PL 
manifold we have constructed is not a differentiable manifold since in that case the transition 
functions should be elevated to a diffeomorphism. That the latter cannot be defined here is 
easy to understand, the intersection of two p-simplices is a (p— 1) simplex, again in the singular 
homology one may introduce the (p — 1) standard simplex which has lower dimensionality, 
therefore the transition function does not have the right differentiability, only (n — 1) partial 
derivatives are defined. 

That we do not have a differentiable manifold is an important limitation when we come to 
consider global objects, such as those associated to vector fields, for example Lie group flows, 
and topological invariants computed not in the Homology (which is trivial using the complex) 
but in the space of forms fi(M; R), such at the Gauss-Bonnet theorem considered in section 
El A look at the space of Whitney forms in 2D (Example I2.1j) will reveal that only the zero 
co-chains are continuous, while the one-cochains are not continuous (See FIO I1.4J) and so are the 
two co-chains. All these are only continuous in the standard simplex A n , i.e locally, by which 
we mean for the given n-degree simplex (Ti( n ). In not having a smooth structure globally, we 
are lead to construct operations locally where there is a smooth structure, as announced the 
interior product will be exact in the sense of locally exact. 

For example, in 2D, suppose the Euler class is discretised as a Whitney two-form, consider 
the integral of the co-chain, say mapping the two-chain e(K) to the Whitney two form e(M) = 
W(e(K)), 

[ e(M) = W e(a), (15) 
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but this is not exact, suppose as in Section|2]we may take the gauge transformation specified by 
P such that in each simplex P~ l 1ZP is flat, then the Euler characteristic is zero. We are missing 
the contribution from distributions over the edges where the 2-form changes discontinuously. So 
one would have for each link a delta function 3 with a weight factor so that e(M) should be 



The factors Kj would be determined by the Gauss-Bonnet theorem, and correspond to the subtle 
play between integration domain and the Euler class which is itself a measure, this interplay 
gives the topological invariant in the continuum theory. In the present case one could imagine 
making a manifold which imitates the present simplicial complex situation in which the field are 
now C°°-functions, in the limit were the vielbein gets close to the simplicial field, we recover a 
delta function as above. Having understood this limitation, namely the lack of relation between 
local and global, we will proceed by means of the overall normalization N for the Euler class 



as will be described in section |3l This seems a round about way to calculate the Euler charac- 
teristic which is x(K) = V — E + F, or using the homology groups which are readily calculable 
in the chain space using the boundary operator d K . 

As seen in FIG. 1, the type of vector field satisfying A introduced is not smooth, for example 
their divergence 



is ill defined over a boundary unless the vector is constant globally because of the hopping of the 
vector field across cells. Also, if the vector field is represented as its dual (in the sense of metric) 
one chain, its divergence is zero in each simplex. It is then clear that the flow of a vector field 
can only be defined locally, it is approximate, but nevertheless the formula for the Lie derivative 
is exact within the cell and compatible with the de Rham cohomology. 

Our aim is to approximate the continuum theory of differential geometry, to that effect we 
construct operations which are locally exact in the sense explained above. On the other hand 
the de Rham co-homology is maintained in the construction, which means that the co-boundary 
operator d K which acts as a linear operator on the chain space, is compatible with the contrac- 
tion for each simplex. In Table ^ we take schemes that may be considered with their list of 
operations and we give on the right hand side the properties satisfied by them as we will see 
shortly. This introduction concludes with the list of the new basic operations and where they 
are defined in the text: 

x, Cartesian product, Eq. liUl 

W, Extended Whitney map, Eq. I41[ 

c/?(°) , Map to function coefficient of a form, Eq. 1421 

D K , Extension of d K to the product space, Eq. EH1 

ir E , Discrete Hodge star on the product space, Eq. 

/\, Discrete wedge on the product space, Eq. |§21 

, Curved metric Hodge star, Eq. I15()[ 
p, Morphism from non-commutative form to usual differential form, Eq. 11921 
T, Truncation map, Eq. 12041 




(16) 




(17) 



V.A 



(18) 



3 The use of delta function for the Euler class is an example of a current |15| 
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Ring of operators 


Properties 


{d K , i x , *, A} 


A, C, P 


{d K , i x , *, A^} 
{d K , i X ,*,/\} 


C 


{d K,L A K,L^K,L } 


A, C 



Table 1. The various schemes and their properties 

2 The scheme for discrete differential geometry 
2.1 Topology 

The mappings between chains and co-chains, W and A 

Consider a hypercubic complex K, for the moment we will work with a square complex. A 
p-chain is a linear combination of hypercubes, 

a (p) = ?.\(p) a i,(p)> ( 19 ) 

i 

the associated co-chain = W(a^) is represented as a differential form, it has well-defined 
value in the standard simplex associated to each n-hypercube of which cr^u,) is a face. One 
associates to the standard simplex the Cartesian coordinate system 

K( P )l = [0,l] p . (20) 

To construct the co-chain as a differential form, the Whitney map and the de Rham co-homology 
inner product are used, the following can be seen as part of the definition of the map W, 



fo,(p)jW( j,(9))) = / %(e)) = ^V ( 21 ) 

The Whitney map is linear, and since d K is the co-boundary operator on the chain space, then 

dW = Wd K . (22) 

We note that the operation d is exact in the interior of every hypercube, this follows directly 
from Eq. 1221 This statement may be expressed as a commutative diagram illustrating the fact 
that W is a functor from the space of chains to the space of co-chains, mapping the morphism 
d K to d. We will refer to this property of exactness as P. 

Example 2.1. Whitney forms for the square complex 

Whitney forms for a 2D square complex In each square choose a vertex as the origin, 
and choose an orientation, then introduce the coordinates (x,y) G [0, l] 2 . Given a square, the 
chain [0123] let [0] be the origin, [1] with coordinates (1,0), [2] with coordinates (1,1) and [3] as 
(0,1). Then, the Whitney zero forms are: 

W([0]) = (1 - x)(l - y), W([l]) = x(l - y), W{[2\) = xy, W([3]) = (1 - x)y, (23) 
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the Whitney one forms: 

W([01)) = (1 - y)dx, W{[12]) = xdy, W([32]) = ydx, W([03]) = (1 - x)dy, (24) 
and the only Whitney two form is 

W([0123]) = dx Ady. (25) 

This is the complete list of co-chains as differential forms, in the chart with coordinates (x, y) 
corresponding to the square [0123], see FIG.\Q 

When expressing the Whitney form in coordinates, it will be convenient to indicate in which 
square it is defined, we will do that by putting the unit function in front of the expression like 

W([32}) = l [0123] ydx + 1 [324 5](1 - y)dx, (26) 

indicating the two squares in FIG.\^ (top degree simplices in which W{[Q1\) has non-zero value). 
This last expression being the global definition ofW{[0\}). 

The de Rham map A, is the left inverse of the Whitney map {AW = 1), it is defined over the 
space of forms, like integration itself, the co-chain space being seen as a subspace of the space 
of forms. For wW G 0^(M; E), define 

A(>)) = ^> >(p) f W W. (27) 
i •'Kim I 

It is immediate that Stokes' theorem may be expressed as 

Ad = d K A. (28) 
The property of exactness is lost as soon as one introduces the discrete wedge product A K . 

The discrete wedge product 

The various mappings are in place now to introduce the discrete wedge product A , it is com- 
binatorial in the sense that it involves the continuum wedge product: 

°i,(p) aK a (Ui) = N (P> l) A (W(a ii(p) ) A W(a Uq) )), (29) 

where the numerical factor N(p, q) is a normalization 4 . It follows immediately that the operation 
satisfies the Leibniz rule, expressed as 

d K (o- ij ( p ) A K <Tj,( q )) = d K cr it ( p ) A (Tj^ q ) + '(^.(p)) d o- j; ( q y (30) 

and is easily proved. The discrete wedge product for the hypercubic complex is better behaved 
than the simplicial one; both are strictly speaking non-associative, but the hypercubic one is 
almost associative in the sense that one may permute the three chains (?i (p), a j,(q)-> a k,(r)i taking 
care of the signs in such a way that the wedge is associative. So we may write symbolically 

a i,(p) aK a j,(q) aK a k,(r) = max(-l) :p o- i) (p) A K ajfq) A K a k r r y (31) 

where V specifies the sign resulting from permuting the chains. The reason for the almost 
associative nature is that in some cases, the result might be zero, hence the reshuffling. 



4 For the square complex it is AT(0,0) = JV(0, 1) = JV(1,0) = 2,iV(l,l) = N(2, 0) = AT(0,2) = 4. Note that 
N(p, q) = N(p + q), so that the normalization is well-defined. 
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Example 2.2. Discrete wedge product 

Label the vertices of a 2D hypercubic complex by a pair of integers i and j, one as integer X 
coordinate and Y coordinates (not to be confused with standard simplex coordinates). We find 
that 

WJ)] ^ K i(i + U)D A* + = 0, (32) 

while 

[(i,j)] A* ([(< + 1, j)] A K + 1, j)]) = + (33) 

is the desired result, hence the need for reshuffling before inserting the brackets. 

The operation of reshuffling up to sign is an anti-automorphism in the space of forms, it does 
not modify the co-homology. 

Discrete Hodge star * , * L 

The last operation to be discussed in this introduction is the Hodge star operator. The procedure 
is to introduce a dual complex L, and the subdivided complex B which consists of all hypercubes 
arising from superimposing L onto K. The Hodge star operator is defined as a map between 
the two complexes (see for example ^5] for more details), 

* K : C (p) (K; R) ^ C [n _ p) (L; R) (34) 



, L 



C( P )(L; R) — > C(„_ p )(-fC; R), (35) 



after introducing a Whitney map for L as well as for K, we introduce the scalar product for the 
Euclidean manifold, defining the norm, 



K(P)>* °i,( g ))= l M W K(P)) AW ( 36 ) 



where W BK and W BL are the expression for the Whitney map on the finer space B. One finds 
that the adjoint 

S K = (-lf(P+ 1 )+ 1 * L d L * K , (37) 

is simply the boundary operator d K , up to sign (see for example |7j). The Hodge star thus 
defined satisfies the defining properties 

5 K = * L d L * K , (38) 

* K * L = (-l) np+ \ (39) 

and respectively interchanging L and K suffixes. Note that the discrete Hodge star is not exact. 
Looking back at Example 12. H we see that W([0]) is as given in the list W([0]) = (1 — y)(l — x), 
while *H / ([0]) = (1 — x)(l — y)dx A dy is not in the list, we discussed this fact already in the 
introductory part. 

2.2 Geometry 

A product on chains is now introduced, it is motivated by the map W used in the discrete 
interior product defined subsequently. The product x will provide us with an exact discrete 
interior product, and we will use it in a way that respects the co-homology. 
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Definition 2.1. Product x and extended Whitney map W 

A product denoted x is introduced on the space of chains \12]l . It is a Cartesian product, the 
ordered pairing 

a i,ip) x a iM) = ( cr i,(p)> <7 .?,(?))- ( 40 ) 
A map 

W : C (p) (K; M) x C (q) (K; E) — ► Q {q) (M; R) (41) 
is introduced, defined as 

W{a it{p) x a m ) = <p<®(W{a ii(p) )) A W{a m ) (42) 

where ip(°) maps a Whitney form to its function coefficient in the standard simplex coordinates. 

Example 2.3. Generic square lattice. 

For a 2D square complex K, see FIG. [H the chain corresponding to a given vertex, say vertex 
[2] which is the chain o"2,(o); gives the form W(o"2,(o))> & ^ as non-zero value over the four squares 
bounding it. 

The edge [32], parallel to the x-axis in the embedding M? gives W([32]) with non-zero value 
over the two squares of which it is a bounding edge. Finally, the square [0123] gives a differential 
form which has non-zero value in that square. 

This limitation of the non-zero support, is one feature of the product x, under the map W, 
we note that the two co-chains must have some overlap where they are both non-zero valued. 
Such considerations will be of crucial importance later, we will have to check that the pair of 
chains appearing in a product are such that the resulting form arising under application of W is 
indeed differentiable, e.g if one considers two edges which only have non-zero value over edges, 
not squares, then the differentiability does no longer holds. 

Lemma 2.1. Generic constant forms 

Introduce the generic constant forms, which can be expressed as a sum of Whitney forms, 

°° = (43) 

(id) 

o 1 = ^[(iJ)(i + lJ)],o*±^[(iJ){i,j + l)], (44) 

(id) (id) 

° 12 = ^2[(i,j)(i + l,j)(i + l,j + m,j + l)}. (45) 
{id) 

They satisfy W(a°) = 1, W{a v ) = dx, W(a 2 ) = dy and W(a 12 ) = dx A dy. 

These forms are globally constant in the sense that within each top degree simplex they are 
constant. 

We note the relation of the extended Whitney map W to the original Whitney map W . The 
following lemma follows immediately 
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Lemma 2.2. Extended Whitney map W 

1. If the chain at p \ is such that W(at p \) is a generic constant form, that is one of the preceding 
equations, then 

W{o- (p) x a iy[q) ) = W(a ii{q) ), (46) 
for any cr i (p) . 

2. The map W acts as the Whitney map W on the chain space 

W{a it(p) ) = W(a it{p) ). (47) 

Here, vector fields live in a Euclidean space. So the dual one form X b to a vector field X 
with respect to the Euclidean metric is represented as a one-chain: 

X K = A(X b ). (48) 

The introduction of the discrete vector field is straightforward by the very fact that topology has 
imposed on us a Euclidean metric. Below, we will see how curved metric data is incorporated. 
Using the chain space to represent the vector field as opposed to the dual complex is we will 
find natural. 

Before introducing the interior product, we need one more map, 
Definition 2.2. Discrete contraction on chains (hypercubic complex) 

Given an edge Cj,(i) and a p-hypercube <Tj,( p ) both faces of a top degree simplex cr k ^ n ^, let 

H a i,{i)i a i,(p)) = i a k,{ P -i) I a k,( P -i) < Ofc.HiOfc.Cp-i) ^ K(T j,( P ) = °) °fc,(p-i) AK(J i,(i) 7^°}) (49) 
then, 

^ l( i)( ff j,(p)) = J2 a k,(p-^y ( 50 ) 

A 

If one draws the hypercube and an edge 0^(1) such that < Then i] a . {1) (o"j,( p )) 

is simply the sum of all (p — 1) faces of orthogonal to crj,(i). 

Definition 2.3. Discrete interior product 

Given a one chain X K , and a p-chain a^, the discrete interior product is defined as 

i : C (1) (K; R) x C (p) (K; M) — , C (1) (K; R) x C ip „ 1} (K; R), (51) 

where 

ixKO(p) = X K x r] xK (a (p) ). (52) 
Example 2.4. Interior product 

Start with a ID example. Consider the edge [01] on the line. Then, the Whitney elements 

are 

{W([0}) = 1 - x, W([l}) = x, W([01}) = dx}, (53) 
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and 



i [01] [01] = [01]x([0] + [l]). (54) 

That is, we took the sum of the two vertices that remain when we remove the edge [01]. Also, 
the construction is guided by the embedding in the continuum: 

W(i [01] [01]) = ^°HW([01})) A W([0] + [1]) = 1. (55) 



The one chain [01] plays the role of a function coefficient so, in this case it is 1 and hence the 
result. If you contract [0] or [1] with the edge [01], you get zero because the edge [01] is not a 
face of either vertices. 

Consider the 2D example of the square (FIG.\^j. This time there are many combinations of 
two chains which lead to a non-zero contraction, for example the contraction. Let us see three 
examples of ix 

W(i [32] [0l23\) = W({32] x ([03] + [12])) = <p<®(ydx) Ady = l [om ydy, (56) 
W(i m [32}) = W([32] x ([3] + [2])) = l [0123] y 2 + 1 [3245] (1 - yf (57) 
W(i [32] [01}) = W([32] x ([0] + [1])) = l [0 i23](l-y)y- (58) 

The second line shows that the form obtained after contraction has non-zero value only in the 
two squares [0123] and [3245]. 

The first calculation shows how when contracting a 2-chain, we get a linear combination of 
the edges perpendicular to the edge which specifies the vector field. This generalizes obviously to 
higher degree chains: we get a linear combination of the (p — l)-simplices which are orthogonal 
to the vector field (an edge). The third line shows that although [32] is not a face of [01], the 
contraction still applies because they are parallel within the same cell of highest degree which 
plays the role of the local open set. 

Application of d K after contraction is prescribed as follows, the extension of d K is defined to 
be D K and acts as 

D K (i X KO- lAp) ) = 2X K x d K ( Vx Ka it(p) ). (59) 
Lemma 2.3. Compatibility of ix with the de Rham co-homology 



The discrete interior product is compatible with the de Rham cohomology in the sense that: 
WD K i X K =dWi X K, (60) 



To prove the lemma, we inspect the list of Whitney elements for the square Eq. (|23|) any 
vertex [i] has within a given top dimensional simplex A n of which it is a face the following 
expression: 

w([i])=n k f^f' l) \ (6i) 

where each factor is either 

f j = ff = (l-x j ), (62) 



or 



fj = *y (63) 
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Assuming that there are p factors of the type / + and (n — r) — (p + 1) factors of the / type 
we get 

w(a (r) ) = n? =1 /; ( . r) n?- +1 /+. r) J\ d x ^, (64) 

k=l 

Then, starting with the left hand side of Eq. EOl we get 

W(2X K x D K V (a)) = d(^((W(X K ))) A W{ VxK (a))), (65) 

By moving the derivatives -t^j across to the left, we pick up all the terms since the vector has 
at least every /j factor that a has. 

Before we proceed we need to make an additional prescription for the wedge product. 

Definition 2.4. Discrete wedge and Hodge star on the product space 
It is defined as follows 

ixo~i(p) A K iyo-j^g) =X ®Y x {r] X ^i{p) A x Vya^) (66) 

where ® is the tensor product. The above formula generalizes easily to the case where there are 
many terms contracted and then wedged. 

A similar property applies to the discrete Hodge star: 



K ' a i(p) X a j(q)) = a i(p) X * a j(q)- ( 67 ) 



Theorem 2.1. The discrete interior product is exact 

W{i X K<J (p) ) = i W (xK)W{a {p) ). 
Furthermore, it satisfies 

(i X K) 2 = 0, (69) 
and it is an anti- derivation 

i X K i a i(p) AK cr j(r)) = ( i X KCr i(p)) ^ K a j(r) + ( _1 ) P<T i(p) A ^ i X Ka j{r)- ( 70 ) 

At this point, we have established the algebraic tools which allow to construct the Lie deriva- 
tive which we take to be 

L x k = i xK D K + D K i X K. (71) 

Furthermore, the Lie derivative is exact, i.e 

WL X K(J i){r) = L wxK W{a i){r) ). (72) 

This immediately implies that one side of the Jacobi identities can be evaluated, namely (by 
composition of the exact mappings) 

[L X , L Y ]. (73) 

To define the Lie bracket, well, it should be exact in order to match Eq. 1731 We note that the 
Jacobi identities for the Lie 

L[ X ,Y] = l L x, L Y ] (74) 
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hold if and only if the Lie bracket is exact, that is if it closes on the chain space 

W([X K , Y K ] K ) = [W(X K ), W{Y K )\. (75) 

However, the bracket may be expressed in terms of the product space. The following can be 
seen as a definition of the bracket: 

{[X K , Y K ] K ) = Z K , (76) 

where Z K is such that 

WZ K = [W(X K ), W(Y K )]. (77) 

Note that the introduction of the product x , enforces the role of the local open set associated 
to the top dimensional simplex, so the various terms arising from the evaluation of Eq. 1741 have 
interpretation in each and every open set associated to a square. 

To think about flows, ax the flow of X. It is well-defined per-cell. However, the vector 
changes abruptly across cells, if X has an orthogonal component to a boundary, it changes 
discontinuously across that boundary. 



2.3 Cohomology 

Before discussing the de Rham co-homology in the context of the product space, we must 
understand how this question arises in the first place. In the picture of co-chains as differential 
forms, as provided by the Whitney W and de Rham A maps, the co-chain have global definition. 
A priori, the product space may contradict this crucial feature since the product of two chains 
a i, (p) x (q) has non-zero value where both cochains have overlapping non-zero value. 

Lemma 2.4. Co-homology and differentiability 

Given a pair of chains (<Tj( p ), °~j,{q)), their product 0i (p) X °~j,(q) is to be acceptable if the 
associated differential form 

W(a i(r) x a j>(p) ). (78) 

is differentiable globally (i.e in each hypercubes). 

i) In such case, the co-homology is not violated by the introduction of the product x. 

ii) Given a one chain and a p-chain, the interior product of the p-chain by the one chain 
gives rise to an acceptable product. 

The statement i) stresses the importance of differentiability and ii) was established in the 
proof of Lemma 12. ^1 



2.4 Alternative definition of the discrete analogues of * and A in the product 

space 

Let us now follow the logic of the introduction of the product x , which was required by exactness. 
The natural question is of its applicability to the other operations in the theory. In fact, we now 
find that if we use x to define a new Hodge star operator, it follows that the operation is exact, 
and there is no need for a dual complex and the Poincare duality holds. 

Definition 2.5. Hodge operation 



The alternative definition of the discrete Hodge star reads 

* SfT i,(p) = a j,(p) x ad(p) 
where d{p) is one of (0, 1, 2, 12) referring to the generic constant forms. 



(79) 
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This definition immediately leads to 
Lemma 2.5. Exactness of the Hodge star 

The discrete Hodge star is exact 

Wir E = *W. (80) 
This statement can be proved by constructing the generic examples, as follows. 
Example 2.5. Discrete Hodge star j{ E 

The simplest example is provided by the constant forms: 

★En o 12 12 /oi\ 
a = a X a = a , (ol) 

The identification is obvious since <f^W(a°) = (p^(l) = 1. Thus, 

W(a° x a 12 ) = 1 A dx A dy = *W(a°) (82) 

so the result is exact. Let us consider the more general case of one edge, i.e not a constant form. 
Now consider 

+ E [(i,j)(i + = [(i,j)(i + 1, j)] x a 2 (83) 

To evaluate this, we need to consider the two squares bounding the edge [(i,j)(i + 1, j)], the 
upper one is denoted A, and the lower one is denoted B. Since we use local coordinates (x,y) 
regardless of the cell used, we introduce the unit function 1a for each square. We wish to show 
that Eq. 12731 is exact, i. e that 

W{+ E [{i,j)(i + = *W([(i,j)(i + l,j)]). (84) 

So, 

W(+ E [{i,j)(i + l,j)]) = <p(°\(l-y)dx)l A Adyl K + <p ( - \ydx)l B Adyl K = (l-y)dyl A + ydyl B , (85) 
which is correct in both squares A and B and is zero everywhere else. 

The proof of the lemma, consists in reiterating the example for zero chains and two chains. 
Lemma 2.6. D K -Cohomology on C^(K; M) 

The chain space is such that under the new Hodge star, 

W+ E C {2) (K; R) < b(°\K; 1 
W+ E C {1) (K;R) < B W {K;M 
W+ E C {0) {K;R) < C^(K;M 

This for the hypercubic complex. 

Clearly this is a limitation if one is to use the chain space rather than the product space of 
chains. Also, consider the property 

S K = (_1)«(P+1)+1* 



(86) 
(87) 

B {2 \K; E). 
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One finds easily that 

S K tr im = (-l)+ E d K + E a m = (90) 

for any 0^(2) • The reason is simple, the Whitney form Wiau^)) is constant in the region 0^(2) 
and therefore since the Hodge star is exact, d K ir E cr^^2) = 0- We are then lead to represent 
fields themselves, in the product space, so not just two chains which are not sufficient as a basis 
for the co-homology with operation 5 defined as above. In general we will consider fields such 
as 

^ K = Yl ^ a h(p) x /%,(?)• ( 91 ) 

Let us now return briefly to the interior product. Because of the exactness of the new Hodge 
star ^ E we can now derive the discrete version of the interior product using Eq (|11|). which was 
approximate in the original construction. To that end, we need to introduce a wedge product 
using the product x . 

Definition 2.6. Discrete wedge f\ 

Given two chains a and p. A wedge product is introduced 

a i,(p) / \Pj,(i) = a hir) ® Pj,(p) x a (p+q)- ( 92 ) 
where is the appropriate constant form of degree (p + q). 

This discrete wedge is numerically exact, by construction, it satisfies 

W(o-i, (r) /\ p m ) = W(cr iM ) A W(a m ). (93) 

We may now turn to the alternative definition of the interior product Eq. f|l 1|) , the approximate 
discrete wedge was used, together with a definition of vector fields in the dual complex. The 
formula is numerically approximate. In the present approach we see that we may follow (101 lllj 
in using the formula Eq. (|11|). but with the exact operations now defined in the product space, 
namely, 

Lemma 2.7. Relation to other formula for the interior product. 

The discrete interior product may be expressed as 

ix*o- ip) = {-l)^-v)+ E X K /\ * £ a (p) . (94) 
It is exact. 

The cohomology is correct by virtue of the exactness of d K applied after the exact operation 
of* s . 

Before turning to the simplicial formulation, we consider briefly the Clifford product as it 
plays a role in the discussion of spinors. The product is 

w v P = ^(-^^few) • • • V) A ia m ■ ■ ■ ia^P, (95) 
p 

this mapping involves both i v and the wedge product. We have defined both maps. Recall that 
the discrete wedge f\ K is approximate so will be the Clifford product |28| . We however found 
some exactness, the discrete version of the Clifford product 

a i,(p) V = X^ -1 ^^!) • • • »a0*) ff i,(p)) AK $mi) • • • ^P^Uq) ( 96 ) 
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is exact if at least one of 0"j ( p ) and <Tj^ q ) is a constant form, i.e is one of cr°, a 1 , a 2 , a 12 . In that 
case it satisfies 

° d v °m = MW(a d ) V W(a m )) (97) 
and is exact. 

2.5 Simplicial formulation 

We now consider the simplicial case. It enables one to consider manifolds such as the sphere, 
but also sheds some light on the method. Now, given one such triangle 0"j ( 2 ), we introduce the 
simplicial version of the Whitney forms in the standard simplex coordinates: 

W([0}) = l-x-y, W([l}) = x, W{[2]) = y. (98) 

A subtle point is that these coordinates do not match the local coordinates of the embedding R N 
up to an isometry. Since differential forms are coordinate independent objects this assignment 
is fine. Moreover, the contracted chain term orthogonal to the vector may be defined in terms 
of the orthogonality in the embedding R N in order to keep geometrical data intact. 
The one forms are: 

W([01]) = (1 - y)dx + xdy, W([02]) = (1 - x)dy + ydx, W([12]) = xdy - ydx. (99) 
The only two-form is 

WQ012]) = dxAdy. (100) 

To construct a valid interior product for the simplicial complex, we need to introduce some 
notation, as there is no direct correspondence between form components (i.e a dx or dy term) 
and co-chain in contrast to the hyper-cubic setting. The first step is to note that the co-chain 
C^ P \K; R) contains the globally constant forms, these are easily constructed, by inspection, 

1 = W([0] + [1] + [2]), (101) 

dx = W([01] - [12]), (102) 

dy = W([02] + [12]), (103) 

dxAdy = W([012]). (104) 

We write the chains on the right-hand side of the equation above as follows 

Definition 2.7. Simplicial notation 

The linear combination of chains which maps to the constant form 

dx J = dx jl A ... A dx jr (105) 

under the Whitney map, is denoted <r/ % . The index J stands for an ordered list of 1^ coordinate 
indices. For example J = {1; 2} gives W(o~^) = dx 1 A dx 2 = dx A dy. 

A co-chain admits various presentations in the product space of chains as for the hypercubic 
version, 

W(a {p) ) = W(a' J xa( p) ). (106) 

the decomposition on the right hand side has the advantage of separating the various form 
components which appear in the simplicial case. 
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Lemma 2.8. Simplicial interior product 
A chain is decomposed as follows 

a {p) = a j x <j( p) . (107) 
The simplicial operation, where X is a one- chain, is 

i x o- {p) = Xfj, ® afu x afp^ (108) 

where a^j is the chain giving the correct coefficient for the W(o"L_^) form component so that 
the discrete interior product ix, is exact: 

Wi x =i W {X)W. (109) 

and it is compatible with the other operations. 



In the formula for the discrete interior product, X^ is the //-component of the chain corre- 
sponding to X, as we have seen, this component can be extracted. Next, the chain a^j appears 
J 

(p-i) 



in the sum a^j x a, where J is the summed index. The chain a^j is such that 



i_|_^(^J x a( p) ) = Wia^j x af^). (110) 
To establish the validity of the formula, we note that 

Lemma 2.9. The space of extended Whitney forms is invariant under contraction by a globally 
constant vector field. 

To show the lemma, we may take the more general case of any dimensions, say N. Then the 
general formula for the Whitney form is usually written using /Xj = W([i]) as a notation, 

W([h . . . i r }) = - { '-IHdfHi A • ■ ■ dfH k A ... A dm,. (Ill) 

k 

Taking the following expression for the Whitney forms, with X{ coordinates in the simplex, 

Hi = Xl , Vi + 0, (112) 

mo = i - y~]xj. (H3) 

i 

With this list, it is found by inspection that the coefficient function of a form is expressible within 
the given simplex as a linear combination of Whitney zero forms, and the form components in 
Eq. 11111 are a linear combination of constant forms. 



3 Further prescription for metric 
3.1 The picture 

Consider the manifold constructed using the standard simplex as the local chart, we find that 
this manifold 6 E has only a local flat metric 5{j corresponding to normal coordinates in curved 
geometry, but if the surface has curvature, the metric is not globally defined in a smooth way, 
it changes discontinuously from simplex to simplex, which are global objects. In the Regge 
approach, the vielbein {e^ a \ a = 1, 2} is associated to edges and the change in the metric 



9 = <We (a) ® e (6) 



(114) 
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from a simplex 0^(2) to a neighbouring simplex cr i+ i^ 2 ) is defined by a matrix, which gives a 
finite rotation of the vector, 

e (a) K (2) ) = A^)((T i+1)(2) ), (115) 

and hence the metric does not vary smoothly, i.e dg is not defined globally, or we would say 
here, that it cannot be represented smoothly as a Whitney form. 

We propose a solution by defining the vielbein smoothly per top dimensionality simplex, 
which means in turns that it is smoothly varying, i.e that within a triangle the metric is not 
locally flat as in the Regge approach. To that end, the product space introduced as a mean to 
represent the discrete interior product will prove crucial. 

Let us start with the manifold structure we are given by the standard simplex construction. 
Given the complex K, we have a piecewise linear manifold with local chart given by the standard 
simplex. The simplicial complex is a geometrical approximation to the manifold S. The discrete 
metric in the product space of chains could be required to approximate this induced metric (it 
need not be the case), and capture local holonomies and other topological features as we will 
see. 

For consistency, integrals will be performed in this chart, so for example: 

vol E (a)= / * B 1= / dxAdy. (116) 

J\a\ J\a\ 

Of course, this Euclidean volume is not the volume in the new metric which is constructed 
algebraically below, the volume in the new metric is 



vole(a) = 




(117) 



We will favour the standard simplex coordinates, and will express 

9^ = 9 x (x,y)dx + 9 y (x,y)dy, (118) 

in terms of Whitney forms. Also, we will use the new Hodge operation it E for * E and at the 
end of this section we will define the discrete analogue of * as it 6 . 

Since the discrete version of the operations {A,d, *, i v } is well-defined and compatible with 
the de Rham co-homology, there is a strong suggestion to use these to build the vielbein metric. 
Ounce an ansatz is proposed for the metric and the Cartan structure equations are written, 
the rest of the construction follows, since the covariant derivative d^ is then fixed, and so is 
the curvature. In the standard simplex coordinates, an equation normalizes the Euler class and 
hence determines partly the metric. 



3.2 Vielbein 

To construct a generic metric, we start with the flat metric 

ds 2 = dx 2 + dy 2 , (119) 

it may be written immediately in the non-coordinate basis 9^ = dx, 9^ = dy. There are two 
different considerations that can be made, namely a change of either i) metric; or a change of 
ii) frame, i.e a change of coordinates, ii) corresponds to a gauge transformation. 

Consider i), since we have chosen a preferred coordinate system, namely the standard sim- 
plex coordinates, which is fixed, it is natural to consider the new metric g = 9^ <g) 9^ a \ where 
0(a) _ P(i x a^ w hii e keeping the coordinate system fixed ( ii) is the alternative). It is of crucial 
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importance that the constant forms dx a may be represented in the simplicial and hypercubic the- 
ories, i.e A(dx a ) £ C^(K) it is not an approximation. Also note that one is deforming a metric 
locally, but as we have seen, the topological features are global and require a normalization. 

For ii), the change of coordinates x' = Px, leads to the new frame e' = eP^ 1 , so that 6' = P9, 
the connection transforms as 

J = P'ujP'- 1 + P'dP'-\ (120) 

while the curvature two form transforms as 

1Z' = P'1ZP'~ 1 . (121) 



A piecewise flat vielbein 

We take the standard simplex coordinates (x, y) as the preferred ones to express the vielbein. 
Let us consider the simplest implementation of metric. In the standard simplex coordinate 
system leads we have the Cartesian frame 

e = R, 8 y }. (122) 

It is then natural to consider a new frame e' defined via a change of coordinates to x' = Px, 
where P is taken to be piecewise constant, the new frame is related to e by 

(123) 

however since the dual basis 6 satisfies 

<0«,e (6) >=S a b . (124) 

it is immediate that 9' = P9. First let us introduce the new coordinates 

a (1) (l-x-y)+/3 (1) :r + 7 (1) y = x' (125) 
a (2) (l-x-y)+/3 (2) x + 7 (2) y = y'. (126) 

To derive the new frame, in (a/, y') coordinates, we derive P form the above equation and fix 
P0 as the vielbein in the old coordinates (x,y), so that it is a new metric (corresponding to i) 
above). 

Since P is constant in the interior of the standard simplex (ft^a, it follows from the first 
Cartan structure equation that 

dP6 = -luAP6 = 0, (127) 

so the connection uj is piecewise constant. Hence the metric is flat in the interior of the cell. We 
have just proved the following 

Lemma 3.1. Family of flat metrics 



The space of chains gives rise to a family of metrics which are flat in the interior of each 
simplex. 

Consider P to be a constant non-singular matrix piecewise defined. Of course this implies 
that the new frame e' gives a piecewise flat metric, with piecewise constant connection, and 
vanishing curvature two-form. 
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Generic metric 

In order to construct a non-zero valued curvature two-form, we need to construct a non-piecewise 
constant P. To do that in a consistent way, we will use the discrete operations introduced above. 
The generic 2D metric is 

ds 2 = dx 2 + G 2 {x,y)dy 2 . (128) 

The vielbein is 9^' = dx, 9^ 2 ' = G(x,y)dy. and so P 1 = 1, P 2 = G(x,y). Immediately one 
obtains the connection and curvature 

W12 = -d x G(x,y)dy, (129) 
H = (d x ) 2 G(x,y)dx Ady. (130) 

The parameters (x, y) are the standard simplex coordinates. Note that the natural choice of 
the vielbein as a one-chain will not do because two derivatives give zero (because it is a Whitney 
one-form), as in the case of constant P. To move forward, we take the following 

Definition 3.1. Ansatz for the vielbein 

Let the discrete vielbein be 

6® = A{dx) (131) 
0(2) = i x + E i Y A{dy). (132) 

The two vectors X and Y are arbitrary one chains, up to giving Euclidean signature, and 
parametrize a family of metrics with line element 

ds 2 = A(dx) ® A(dx) + i x + E i Y A(dy) <g> i x * B iyA(ciy). (133) 

To obtain an orthogonal metric we require that X and Y be such that 

ix+ E iYdy = G K x A(dy), (134) 

for some chain G K to be determined, since dx and dy are orthogonal in the standard simplex 
coordinates. The ansatz for the vielbien should be understood in terms of connection and 
curvature. 

Lemma 3.2. Discretized Cartan structure equations 

The discretized version of the torsionless Cartan structure equations is given by 

d K 9 (a) +i( ,E^ Kir E 9 (b) = Q) (135) 
dV + *V*V = n K . (136) 

The ansatz allows to solve the first equation, the second is a definition of the curvature two- 
form, the metric is normalized by the Gauss-Bonnet theorem. 

To prove the lemma we will consider a generic example below, first we solve the problem in 
the continuum theory by calculating u. Now, choose X and Y such that Eq. (J265|) holds. Next, 

d9^ = {d x X x Y y + X x d x Y y )dx A dy, (137) 
since the second derivative (d x ) 2 gives zero on chains X and Y, we are left with 

W2i = ~{d x X x Y y + X x d x Y y )dx. (138) 
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Next, in 2D the curvature two- form is simply 

■R = duj = 2d x X x d x Y y dx A dy. (139) 
Finally, the volume form is 

/l = 0WA^. (140) 

The curvature two- form 1Z is piecewise constant, it is therefore a two-chain and does not require 
the product space to be represented. 

Up to now we have defined the discrete objects, and solved the continuum version of the first 
Cartan structure equation for the connection. We now introduce the covariant derivative. 

Definition 3.2. Covariant derivative and curvature 

Define the covariant derivative d u , as 

d w K = d K + (141) 
Then the curvature two-form is given by 

Tl K {u K ) = (d wK ) 2 (142) 
mimicing the continuum formula lZ(u>) = d 2 ^ and Eq. (|136|) . 

We note immediately that 
Lemma 3.3. Exactness of d u and 1Z 

The discrete covariant derivative is exact and can be expressed as 

Wd u K =d w ^ K) . (143) 
In turn, the curvature is also exact: 

K{W{uj k )) = W{U K {uj K )). (144) 

The lemma implies that the curvature two form is the true curvature corresponding to the 
connection u . 

To consider parallel transport, 

d^P = 0, (145) 

one may write equations similar to the first Cartan equation. 

To determine X, Y in the ansatz for the vielbein, we find that these are constrained by 
topology, we must impose the Gauss-Bonnet theorem, depending on the number of independent 
variables Y, this is a constraint that does not determine the volume of the surface in the new 
metric. 

To make contact with topology, we must calculate the first Pontryagin class p\{K) and show 
that the Euler class gives the correct Euler characteristic upon integration. Now, 

MK) = ^trTl 2 . (146) 

but we have computed explicitly P, so the curvature is modified to P~ 1 1ZP, while if we con- 
sider the spherical metric ds 2 = (d8) 2 + (sin 6) 2 (d(j)) 2 , we have a new P = P' such that 
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1Z = P'~ l 1Z S phereP' \ and by cyclicity of the trace, we get the same result, but only locally, 
i.e as we discussed, this argument does not hold globally and so we appeal to a normalization 
for the Euler class. We return to this point in the next section where an explicit spherical 
topology example is worked out. 

Let us now turn to a simple example, to illustrate the method in more detail. 

Having introduced a map ir E which corresponds to the Euclidean Hodge star, we are now lead 
to consider the discrete Hodge star for the generic metric if . After application of the curved 
metric Hodge operator one obtains a multiplicative y/g factor, this factor is a determinant and 
determines an orientation like the volume element on a Riemannian manifold 



which has " outer orientation" . 

We have proposed a new discrete Hodge star which extends the former one for a curved 
metric which does not require the introduction of the dual lattice and handles the Hodge star 
in a way that has the correct algebraic properties. We now proceed to describe it in details. We 
start with the volume element 



vol = \fg~dx 1 A ... A dx N , 



(147) 



vol = e (0) A... A0 (3) . 



(148) 



After discretisation, it may be written as 



vol K = W(G K xa 012 ). 



(149) 



In that case the factor G is the discrete analogue of y/g. The discrete Hodge star is then given 
by a map from the product space of chains onto itself. 

Definition 3.3. Discrete Hodge star it 9 



+ 8 = G K 




(150) 



And 




is only formally defined 5 in the following sense: 



(151) 



WiiG*)- 1 ) = (WiG*))- 1 



(152) 



One verifies that indeed 



^(★w 1 ) = ^((Vx* 9 )- 1 ) 

= W(G K )W((G K )- 1 ) 



(153) 
(154) 
(155) 



= 1. 



We also verify that the new Hodge star is related to the vielbein basis, 

Vi = ^a^ 2 ),*V 1) a^ 2 ) = i 



(156) 
(157) 



Now, the Hodge star is primarily used to define the coboundary operator: 



(5 = *-^*. (158) 

5 It cannot be represented exactly as a chain or product chain, due to the functional form of Whitney elements. 
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Which in the present formalism should be discretised as 

8 e = (G^)- 1 x + E dG K x * £ . (159) 

In this context, the discrete exterior derivative is piecewise exact, while the co-boundary operator 
is piecewise approximate for curved metric, since for the flat case, the Euclidean Hodge star ir E 
is exact and in that case (ir E )~ 1 = it E . 

To establish the identity (5 6 ) 2 = 0, we calculate 

(5 e f = (G K Y l x + E dG K x + E {G K y l x + E dG K x * B (160) 
= (161) 

Next, let us turn to the curved metric Laplacian. For a general metric, in the continuum, on 
zero forms for example 

M = ~d vy /gfd v f. (162) 
V 9 

It is specified in the ring of operations we have introduced as 

A K = D K 5 K + 5 K D K , (163) 

where 

S K = + e D K ir d . (164) 

3.3 Explicit construction for a complex with spherical topology 
Choosing a discrete metric parametrization X and Y 

For simplicity we take the boundary of the tetrahedron as a triangulation of S 2 . The faces are 
[012], [213], [320], [013]. We make an ansatz for the vectors X,Y: 

y = £**m (165) 

X = i dy YxA(dy). (166) 

ans use the general formula given above for the discrete vielbein. We also introduce a normal- 
ization for the vielbein 9^ = N6^ a \ This normalization does not appear in the formula for the 
spin connection uj as given in the Cartan structure equation. It will allow us to fix the overall 
volume of the sphere. 

Next, using the explicit form of the Whitney one- forms, we find 

Y y = Y 01 x + Y 02 (l-x)+Y 12 x, (167) 
Y x = Y m {l-y)+Y m y-Y l2 y. (168) 

By direct computation, this leads to 

q{2) = N ( Yoi x + Y 02 (l-x) + Y 12 x) 2 dy (169) 
U2i = -2(Y 01 x + Y 02 (l-x)+Y 12 x)(Y 02 -Y 01 -Y 12 )dy (170) 
n = -2(-Y i + Y Q2 - Y 12 fdx A dy. (171) 

We are computing as in the continuum theory, but with Whitney forms, because the formula 
for the discrete Cartan structure equation and exactness of d K allows such identification. 
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Killing vectors 

The equation is 

L z g = 0. (172) 
That is 

L z 9^ + fl« <g> L z 0« + L z 9^ ® 6^ + 0V> ® L z 6® = 0. (173) 

So let Z = and let us look at the triangle [012], then 

L z 9^ = (-Z 01 + Z 02 -Z 12 )dy, (174) 
L Z 6W = 2Y y Z x (Y 01 -Y 02 + Y 12 )dy + (Z 01 -Z 02 + Z 12 )Y y 2 dx (175) 

Where we follow the notation Z x , Z y as in Eq. (|167j) . The Killing equation are then 

-Z 01 + Z 02 -Z 12 = 2Y y Z x (Y 01 -Y 02 + Y 12 ) (176) 
(Z i - Z 02 + Z l2 )Y y 2 = (177) 

These equation cannot be solved exactly as they stand apart from trivial case. The de Rham A 
can be applied to the components and then solved, otherwise we may minimize the expression 
WL z g. Also note that there are such equations as Eq. I|174|) . (|175|) for each simplex of the 
triangulation. 

Euler characteristic 

To address the topological invariant we consider the simplest case we let Yoi = ^02 = Y12 = Y ■ 
Then, 

fl( 2 ) = (Y + Yxfdy, uj X2 = 2(Y + Yx)Ydy, K = 2Y 2 dx A dy. (178) 

Then, 

pi = = ^2 ( Y2dx A d y) 2 ( 179 ) 

so the Euler class is the "square root" of this last expression, 

e(M) = -^-=Y 2 dx A dy, (180) 

7T\/2 

it is a volume element itself, and satisfies the Gauss-Bonnet theorem 

/ e(M) = X (S 2 ) = 2. (181) 

This determines Y, by integration, where Nt is the number of triangles. So we have the following 
equation from topology, 

E / e ( M ) = = 2- (182) 

where the area is calculated in the standard simplex coordinates as 

f dx f dy = \. (183) 
Jo Jo z 



27 



So we found that 



We may now evaluate the volume of the surface S with the metric {0^}, 
vok(K)=J2[ *«A*« = ^£ 



2 r 287T 2 

(l+x) 2 dxAdy = — -. (185) 

i,w 3Nt 



To limit the discussion to our example of one tetrahedron, Nt = 4 and we find 

vole(K) = (^-) 2 . (186) 

Also local holonomies may be considered, by taking a given triangle [012]. In its interior we 
write, what is simply the Stokes theorem (this is the 2D case, in general the exponential is used 
to write the holonomy) 



K = / u. (187) 

| [012] | J\d[012}\ 

This concludes the discussion of metric. 



4 Relation to the non- commutative geometry 

As established by Woronowicz Jl], a non-commutative differential geometry is associated to a 
graph (see also |24j). it is a deformed calculus with deformation parameter of the order of the 
lattice spacing. The typical graph that determines such calculus is a Hasse diagram, i.e a graph 
with arrows. 

Consider the commutative group lattice C (the vertices of C coincide with the vertices of 
the hypercubic complex K), where the group operation is given by translation by one lattice 
spacing. A Hopf algebra can be constructed in such system (in the context of physics see \25\). 
and the Yang-Baxter equations are found in terms of the discrete wedge product, concretely, as 
we recall below, the non-commutativity is between differentials and functions. To start, Let us 
now compare the functions on the group lattice with the Whitney zero- forms. 

4.1 Functions 

The basis of functions in the non-commutative formulation is {e x , x G £}, where x labels a 
vertex of the group lattice C The algebra of functions is 

e x (y) = 5*y, (188) 
e x e y ± § xy e x (189) 

A trivial property is eX = O n the other hand, the Whitney zero forms, {/Xj, i £ K} are 
such that 

MiO') = <% (190) 

which corresponds to Eq. (|188f) and also satisfy J2i fa = 1- 
The second property Eq. (|189|) may be written as 6 

MA*b1=<%b1 (191) 



6 A vertex in the non-commutative picture is denoted x, in the present formalism it is denoted by i as an index 
and by [i] for a chain. 
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This gives a natural correspondence with the definition on a group lattice Eq. ()189|) . We will 
see shortly that such an analogy is limited to zero forms. Next, we introduce a morphism p, by 
first setting 

p : £^C (0) (K,R) (192) 
p : e x ^ m. (193) 

where in this case x and [i] coincide. Now we define 

p(e x e y ) = 5 xy p(e y ). (194) 
which make p a mapping to the space of forms. 

4.2 Forms and the morphism p 

To develop further the correspondence we need to consider higher forms. At this point the 
non-commutativity between functions and forms appears in the non-commutative theory. We 
specialise to the symmetric lattice [21] for reasons that will become clear below. 
The basis list of one forms consists of 

9» = ^Vde 3 ^, (195) 

X 

0"" = ^Vcfe 3 ^, (196) 

X 

where p is the unit vector in the p direction on the group lattice. This choice is special in the 
sense that it keeps locality, e x and e x=t ^ are nearest neighbors. To define higher degree forms 
one needs to define the wedge product for this non-commutative basis, 

6* A 0" = 0» <g> 9 V - d u+t " (8) 0*. (197) 

This operation leads to the braid matrix (see for example ) . In the sequel we drop the symbol 
(g) for the tensor product of non-commutative forms. 

The non-commutativity property of the product of function and one form is 

e xTVQ±H = Q±H e x_ ( 198 ) 

This algebra of functions and forms is non-commutative so the correspondence with de Rham 
differential forms is apparently lost since the order in which the various terms are presented 
corresponds to different differential forms. To remove this ambiguity, we note that there we can 
choose a preferred presentation of non-commutative forms that is unique 

uj= ^ e x *°de x *i ...de x ^. (199) 

il<...<i r 

We now make a prescription for de Xi p by extending the range of p from functions to forms. Since 
in the non-commutative setting the discrete exterior derivative denoted d satisfies the Leibniz 
rule 

d(e x e y ) = de x e y + e x de y , (200) 
we extend the morphism p to accommodate this property, 

p(d(e x e y )) = d(p(e x )p(e y )) = d(pipj) = dpipj + pidpj. (201) 
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The last equal sign is the continuum Leibniz rule on Whitney zero- forms. So we may now write 
the correspondence for any form, using the functorial property of d, we simply set 

p(%2 e x *o de x ^ . . . de x ^ ) = ^ p jo dp h A ... A dp jr . (202) 

i 3 

So we have found that under the map p, maintaining the Leibniz rule leads to the prescription 

p(de x ) = dp(e x ). (203) 

Note that the form of Eq. 1)202(1 is reminiscent of the simplicial Whitney form expression, Eq. 
piljl but it does not correspond to the hyper-cubic Whitney form. To establish a relation we 
introduce a new map T in the following 

Lemma 4.1. The truncation map T 

For a hyper- cubic complex, the Whitney form may be expressed as: 

%)) =T^{-l) k p Ji dp Jio A ... A dfH k A ... A dp ir (204) 
k 

where T acts on the function coefficient, any factor in the expression Eq. ()204j) is of the form 
(x fl ) K or (1 — x^) , where K is a positive integer. The truncation map is such that each term 
of the form x^ or (1 — x^) shows ounce, i.e any K is set to unity. 

To illustrate the map T, choose a hypercube, and maps its vertices to pi functions, for 
example [01], 

rOMMi-jMMo) = T{{1 - x){l - y)d{{x{l - y)) - x(l - y)d{{\ - x){l - y))) (205) 
= (1 - y)dx, (206) 

which is indeed VF([01]). So the truncation map allows to recover a Whitney form. This allows 
to show the relation between the non-commutative and the commutative settings. The maps 
are 

n(£; E) n(M; R) 

\ 4> \T (207) 
C(K; R) 

where 

<p = Top, (208) 

the point is that the composition of the two maps is well-defined and hence the relation to 
the non-commutative basis can be made. The space 0(£; M) is the non-commutative one. To 
construct the non-commutative space of forms, we start from the co-chain space and invert the 
maps. In 2D, zero forms in the square [0123] are given by 

Mo = (l-x)(l-y) = Tpe°, (209) 

/ii = x(l-y) =Tpe\ (210) 

P2 = xy = Tpe 2 , (211) 

p 3 = {l-x)y = Tpe\ (212) 
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Where T plays no role as no truncation is needed. Next, for one-forms recall the example above. 

W([01]) = (l-y)dx = T(l-y) 2 dx = T(pe°de 1 -e 1 de°), (213) 

W([12}) = xdy = Tx 2 dy = Tp(e 1 de 2 -e 2 de 1 ), (214) 

W([32]) = ydx = Ty 2 dx = Tp(e 3 de 2 -e 2 de 3 ), (215) 

W([03}) = (l-x)dy = T(l-x) 2 dy = Tp(e°de 3 -e 3 de°), (216) 

and the only two-form is 

WQ0123]) = -T(W([01]) A W([03]) + W([32]) A W([12]) + W([01]) A W([32])) (217) 
3 

The important point here is that the truncation map projects a general form spanned by a 
subspace of f2(£; M) onto a subspace of forms which matches co-chains. The truncation map 
is invertible, i.e we can recover the form that is in one-to one relation with a non-commutative 
form (see Eq. (127721 ). 

4.3 Operations related to contraction 

The way we constructed the contraction was to involve a product space of chains. Such details 
are not needed when the contraction only involves constant vector fields, namely the dual "^jr of 
dx^ in the Euclidean embedding space. We write the interior product as vm which is identified 
with i_a_. In that case, given a chain 07„\, we find 

dxf 1 

V"0-(p) =T7<7f(o-(p)) (218) 

which is a linear combination of (p — 1) chains. The prescription for the relation to the non- 
commutative setting (see for example |27j ) is the following 

»7«r"(<7(p)) = P( e t*P*(<7(p))), ( 219 ) 

where p* is the pull-back map, it associates a unique non-commutative differential form to a 
chain. The map of the non-commutative theory is 



m " 7^77, r-„ +t^-7 T/i (220) 



- | _d 
where 

T± fl f([(i,j)] = f([(i,j)±m- (221) 

The displacement operator being inserted to compensate for the non-commutativity involved in 
swapping the function and form before operating the derivative JL It has been shown that the 
contracted form in the non-commutative setting is free of displacement operators. This is a hint 
of the matching we expect in Eq. Q219|) . 

We now move to the discussion of the Clifford product. Consider the square [0123], and local 
coordinate chart (x,y). We have four functions e°,...,e 3 and four one forms de°,...,de 3 . 
Consider the form e°6 x . Then, 

idA^odpi - pidpo) = i 9x (l - y)dx = (1 - y) = p + p 1 . (222) 

find that the functional form given by the right hand side of ()202j) is preserved. 
Consider two one forms f^dx^ 1 and g v dx v , their Clifford product is 

V g v dx v = f^gJw + f^dx" A dx v (223) 
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For example, 

W([01}) V W([01}) = ^^([01]) A ^^([01]) = (1 - y) 2 . (224) 

This is problematic, so let us pull it back to the non-commutative form basis: 

-y)(l-y) = (e° + e 1 )(e° + e 1 ) = (e° + e 1 ) (225) 

We have to check that this truncation of the function is consistent with the associativity property 
of the Clifford product. Let us check the associativity with an example: 

Example 4.1. Associativity of the Clifford product 

p*T*([01] V K [02]) = p*T*([01] A K [02]) (226) 

= {e°de l -e l de°){e°de 2 -e^de 1 ) (227) 

= e°de 1 de 2 . (228) 

p*T*([02] y K [03]) = e°de 2 de 3 . (229) 

T/ien, 

[01] ([02] [03]) = e°de 1 de 2 de 3 , (230) 

([01] V K [02]) V K [03] = e°de 1 de 2 de 3 . (231) 

In general, there are two cases for the multiplication of functions as follows, the first one is 
when we are multiplying two identical functions like (e° -1-e 1 ), in this case a truncation is needed 
i.e T(l — y) 2 = (1 — y), the other instance is when there are different (e 1 + e 2 ) multiplied by 
(e 1 + e 3 ) then there is no truncation i.e 

(e 1 + e 2 )(e 1 + e 3 ) = e 1 (232) 

is the correct result after mapping both sides under p. A special case is when we take the 
Clifford product of two edges which are parallel within the cell. For example [01] and [32]. In 
the original picture, 

T*([01] \/ K [32}) = (l-y)y, (233) 



p*((l - y)y) = ((e° + e 1 )(e 3 + e 2 )) = 0. (234) 
Thus this destroy associativity? 

To settle this point we take the edge connecting the two parallel ones, [03] and consider the 
associativity. 

p*T*([01] y K ([32] y K [03])) = (e° + e 1 )(e 3 ) = (235) 
which is equal to 

/T*(([01] V K [32]) V K [03]) = ((e° + e 1 )(e 3 + e 2 ))(e° + e 3 )) = (236) 

So the associativity is kept. The examples show that the correspondence between the Clifford 
product and its non-commutative version is as follows: 
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Lemma 4.2. The non- commutative Clifford product V is related to the discrete product V K by 
application of the truncation map T, 

p*T*(a V K rj) = p*((TV) V (T*7?)) (237) 

the right hand side is an associative product, therefore so is the left hand side. 

What this means for the relation to non-commutativity is that the operations of contraction 
and Clifford product defined here admit a well defined projection to the non-commutative op- 
erations on a group lattice, however the correspondence involves a truncation of the function 
space which violates the algebraic properties of the operations. 

This partial failure of the correspondence for the latest operations is conceptually no surprise. 
There is no satisfactory definition of these operations in the non-commutative picture. This 
applies also to the Hodge star. That is using the non-commutative geometry as a substitute for 
the ordinary differential calculus does not work. 

5 Discussion 

In order to test the usefulness of the method in concrete applications, we propose three concrete 
applications: 

• Topological field theory: calculation of the Reidemeister torsion in a combinatorial way. 

• Simplicial gravity: comparing with 2D random surface models for gravity, for Nambu-Goto 
and for actions including matter fields. 

• Electromagnetic theory: discretising Maxwell's equations. 

The first application is concerned with the calculation of the Redeimeister torsion |33j by ex- 
plicitly constructing a simplicial regularization of the Abelian Chern-Simons theory. In ^B] 
the geometric discretisation scheme using the dual complex based discrete Hodge star * , * L 
recalled above was used. It was shown in this context that 



-J— det(dido) det^di)- 1 det(d 3 d 2 ). (238) 

where Nq is the number of vertices and N3 the number of three simplices of the simplicial 
triangulation of the lens space or of S 13 as discussed in the above paper. The overall factor 

jy"jv is then the metric dependent factor of the expression for the combinatorial torsion. 
This factor can be either derived from the original expression of Ray-Singer, or can be deduced 
numerically when showing the subdivision invariance of the torsion. 

In the present formulation, we expect the following functional form for the torsion 

T' = ^Fdet(«5f D«)\z* c{K . R) det(5f Df)^! c(Ar . R) det(6« D«)\^ C{K . R) . (239) 

We have replaced the discrete operators of the geometric discretisation by their new version on 
the product space of chains, and we have added a factor F to be specified, also note that each 
determinant is restricted to a certain subspace of the product space of chains. To relate the two 
quantities, we relate the operators. Let us start with the first determinant term det(9i<io) of T 
and the corresponding one for T . 

The key is the relation between the chain space C(K;M) and the product space of chains. 
The latter is not a chain space, so we need to related it to the chain space C(K;M). First we 
consider the mapping 

$ = AW : C {0) {K;R) x C {p) (K;R) — ► C {p) {K;R). (240) 
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By inspection, the space ker $ is not empty. An example is a product chain such as [0] x [12]. 
Next we introduce a pull-back map which gives a linear combination of product chains, 
this pull-back should be carefully chosen so that the equations below indeed hold (an example 
follows). The space <&*C(K;M.) is the topological representative of the chain complex in the 
product space. It can be shown that for any chain f^p), 

3>*(d« +1 d«a i>(p) ) = M^ +l D^^{a h{p) ). (241) 
The matrix A corresponds to fixed degree dependent numerical factors. It follows that 

det(d^ +1 d^) = det^Ap+K^L^V) (242) 



/ 



= det(A) det^j-D*) (243) 

where det' is the determinant in the ^*C^(K; R) basis denoted as det(.)\$*c(K;R) m Eq. 12391 
Repeated application of the argument suggests that T = T' if 

F = IvHvT det ( A 3) det ( A 2) det(Ai). (244) 

Since this argument is somewhat abstract, we proceed with a calculation that illustrates the 
above calculations. Consider a single tetrahedron as a discretisation of S 3 . We have shown 
that the operations on the product space are exact. To evaluate Sd, we take one element of 
the chain basis [12] G C^(K;M). In the product space we take Ai[l] x [12] and A 2 [2] x [12] G 
C(q\(K;M) x Cm (_RT; R) . First we carry out the calculation in CVy (_RT; R) , it reads, 

d K [12] = [312] + [012], (245) 
d K d K [12] = [12] - [32] + [31] + [12] - [02] + [01]. (246) 

Now let us evaluate the analogue in the product space as 

W([l] X [12]) = x(xdy-ydx) (247) 

6W([1] X [12]) = 3xdyAdy (248) 

dSW([l] x [12]) = -3dy = 3W([02] - [12] - [23]). (249) 

Next, 

W{[2] x [12]) = y(xdy-ydx) (250) 

dSW{[l] x [12]) = 3dx = 3W([01] - [12] - [13]). (251) 

By comparison of the calculation in the chain space with that in the product space, we find that 
if Ai = — A2, and adding the two terms, the product chain 

a = A([l] - [2]) x [12] (252) 

is such that 

A(6dW(a)) = 3Xd K d K [12}. (253) 

This factor of three will appear as an entry of the matrix A. What we have seen is how the 
chain space is captured in the product space and how the new operators 5 K and D allow to 
compute the torsion. This needs to be confirmed numerically, furthermore, we may compute 
the determinants above in the entire product space to find the overall geometrical factor. 

The second application is a test of the metric aspects of the theory and of the ability to 
express actions with the method. The theory of 2D gravity (string world sheet trajectories) 
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is well understood, and a combinatorial solution of the problem has been constructed (see e.g 
[2HD- It therefore provides a good test for discretising gravity before going to higher dimensions. 
To formulate a valid analogue of the path integral formalism in this theory, we first define the 
ensemble of all paths as the set of all equivalence classes (under diffeomorphism) of surfaces X: 
S 1 I — ► K d with boundaries 71, 72 and we choose the action of a surface equal to its area: 

S NG (X) = J dA(X)= j d 2 £Vdet h (254) 

dA(X) denotes the area element of the surface in and h is the induced metric given by 
h a b = 75p"75p"> where (^ 1 ,^ 2 ) are local coordinates on the manifold S 1 ^) I. Using the Regge 
calculus, the Euler characteristic can be readily computed. In our approach, the metric is not 
the induced metric by embedding, as will discuss below. The partition function is defined by: 

<3(7i,72)=/ D[X]exp{-S NG{x) ) (255) 
J[X] 

Where the integral is over the space of equivalence classes of surfaces. Consider a piecewise 
linear surface with action with chemical potential fj,, where Nt is the number of triangles, the 
action is 

S = pN t . (256) 
The partition function is: 

Z(») = Y,-^T^P(-» N t(T)) (257) 
Ter T 

where Ct is a symmetry factor of the triangulation. We wish to calculate expectation value 
of a physical observable II. To that end, the Monte Carlo algorithm creates a sequence of A 
statistically independent configurations T a with a probability distribution: 

P(Ta) = ^75^ eX P(-^( T ")) ( 258 ) 

and evaluates the expectation as: 

j a=A 

<n>=-^n[T a ] (259) 
0=1 

where A is large. We expect a strict equality as A — > 00 The configurations are generated by a 
Markov chain t which can be regarded as a random walk in configuration space with transition 
function 

i f P \T b ] 1 

T a ,T b \ = Cmin|l,-^M (260) 

The calculation of observables such as the Hausdorff dimension of critical surfaces in the theory, 
is done by a Monte-Carlo simulation, or importance sampling. For a given genus, an algorithm 
generates random surfaces in the ensemble of surfaces with same topology, as determined by a 
probability for an update which is a functional of the action S. The latter will in general be 
expressed in terms of the Regge calculus. 

In the present case one proceeds differently. The metric provided on the surface is not the 
induced metric from embedding. The 2D metric is given by the prescription given in the section 
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covering the metric above. Instead of modifying the surface by update as in the former method 
we propose to update the metric assignment {Yij}. In this case the NG action reads 

S' NG = J 0« A 6 {2 \ (261) 

The discrete version of it can readily be evaluated. It is 

S NG = ^2 f W(6 {1)K ) A W(6 (2)K ). (262) 

To proceed, in each simplex one must assign an origin together with x and y axis. These are 
naturally associated to vertices in the standard simplex construction at work here. So for the 
two simplex 0^(2), We relabel its three vertices as 0(i), x(i) and y{i) respectively. 
The discretized action is 

S = \ Y 0(i)y(i) + J^ Y 0(i)x(i) + ^(Xo(i)x(i)Yx(i)y(i) + Yo(i)y(i)Yo(i)x(i) + Y x(i)y(i)) (263) 
i 

The update in the ensemble of surfaces is now replace by an update in the ensemble of metrics 
parametrized by SYij = ±1 for a fixed underlying surface. We will denote the updated metric as 
Y i ^ Y(5Y). After randomly choosing the edge [ij], the probability for the update is a function 
of the corresponding change in the action. 

At each update, one must make sure that the Euler characteristic still holds. This can be 
done by updating the edges which have one vertex in common with either [i] and [j] in such 
a way that the topological identity still holds. To be precise, when the metric configuration is 
initialized, it is normalized so that the Euler characteristic is verified by integration of the Euler 
class. For the generic 2D metric the Euler class is 

e ( K ) = —m Yl(~ Y 0(i)x(i) + Y 0(i)y(i) + Y x (i) y (i)) 2 a iji2 ) (264) 
the constraint from topology is then 

J W(e(K)) = 2. (265) 

At the first step, an edge is picked, Y^. Next, we find each two-simplex of which it is an edge, 
making the list of two simplices k G / if [ij] < Cfc,(2)- In turn we take all edges which are 
boundary of these simplices apart from [ij] itself. The value of Y is updated for each of these 
edges such that the Euler class is left invariant each such edge is update with the same value 5Y', 
after the two updates Y i-> Y(5Y,5Y'). The invariance of the Euler class means the following 
integral equation 



/ 



W(e(K)(Y)) = J W(e{K)(Y(5Y,5Y'))). (266) 



The update 5Y' is a topological correction to the metric update SY. 

The advantage of the present approach, with its systematic discretisation of differential ge- 
ometry is that actions (at least bare actions) for matter fields can be readily constructed as 



Sm 



= j W(<f> K ) A W(+ e D K <S> K ) (267) 
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where $> K is some physical field and D K the differential operator. In the case of fermionic fields 
one must keep in mind the doubling problem and so one may consider including a Wilson type 
term using the Laplacian A . 

As the third application considered, we are interested in the discretisation of Maxwell's 
equations (see e.g |31j). In the gauge theory formalism the field strength is 

F = dA, (268) 
which in terms of the physical fields E and B is 

F = B + E Adt. (269) 
Then, Maxwell's equations are simply 

dF = (270) 

SF = j. (271) 

In the 4D formalism, the interior product is used to express the boundary conditions which may 
be written as [32] 

in(n A (F 2 - FJ) = (272) 
i„(nA(*F 2 -*iq)) = j s (273) 
where the subscript s refers to surface 

j s = Ps -J S A dt. (274) 

and the vector n = where / is a function such that df is orthogonal to the boundary. 
The present formulation may be used to access such situation in the discrete setting. A full 
implementation is needed to confirm this point. 

On the conceptual side, we note that in the literature a distinction is made between forms 
which are represented by chains and twisted forms represented in the dual complex. This has 
much to do with the Hodge star operator, if one look at the (3 + 1)D formulation, i.e the 
discretisation of the spatial part of the Maxwell's equations, we have the metric dependent part 
of the Maxwell's equations given by 

D = e*E, (275) 

B = fi*H. (276) 

The fields E and H are arguably forms and the fields D and B are twisted forms. Simplicially, 
one argues that twisted forms have an outer orientation, in fact this is just the situation of a 
volume form or density in the continuum. In the continuum, one may introduce a pairing to 
represent such twisted forms as (a, Q) where a is a form and £1 is a volume element |3()| . This 
pairing is interesting here since it can be identified with the pairing x central to this paper, 
so that the latter appears also as a discretisation of the pairing of Burke in the context of the 
electromagnetic theory. 



6 Conclusion 

A discretisation scheme has been proposed in which some key aspects of differential geometry 
relevant to physics are captured. A product space of chains introduced here and the extension of 
the Whitney map for the product space played a central role in the construction. Furthermore, 
the vielbein formulation of metric is consistently discretized in a smooth way and topological 
invariants are used to normalize some parameters in the metric. For comparison we established 
a correspondence, but no identification with the non-commutative geometry of graphs which we 
argued is incomplete as to capturing continuum like metric data. Three applications to physics 
are sketched, and are the focus of ongoing work. 
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Figure 1. The two dimensional plane with a square complex, coordinates (x, y) and a Whitney form 
with function coefficient sketched in the W-ax\s. In FIG. it corresponds to the co-chain VF(2[03] + 
[12] + 3[35] + [24] + 2[67] + 3[78]) 
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Figure 2. Part of a square complex, with the one-chain X corresponding to the vector field 



